Spin-coupled double- quantum- dot behavior inside a single- molecule transistor 
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We report on the observation of Kondo and split Kondo peaks in single-molecule transistors 
containing a single spin transition molecule with a Fe 2+ ion. Coulomb blockade characteristics 
reveal a double quantum dot behavior in a parallel configuration, making our system a molecular 
equivalent to a semiconducting double-quantum-dot system. As the gate voltage is increased the 
charging of the second dot by an additional electron induces a splitting of the Kondo peak. We 
discuss possible origins of this effect including a spin transition into a high-spin state. 



The Kondo effect, originally observed in diluted mag- 
netic alloys is induced by the interaction between a lo- 
cal spin and conduction electrons jlj. In a quantum 
dot(QD), this effect gives rise to a zero-bias conductance 
peak (ZBP) for charge states with a total spin S=l/2 
[2, 3]. When a magnetic field is applied, the degeneracy of 
the spin state is lifted by the Zeeman effect and a splitting 
of the ZBP is observed [4] . In the absence of a magnetic 
filed, a split ZBP corresponding to a bias-induced tran- 
sition from singlet to triplet state can also be observed 
for charge states with an integer spin [l-Q • Another sys- 
tem where the Kondo effect has been investigated is the 
two impurities system originally described by Jones et al. 
0| . For example, split ZBPs have been observed in dou- 
ble QD systems where two dots interact via a conducting 
region or via direct tunneling 0, [To[ . The exchange inter- 
action between the spin of the two dots lifts the degener- 
acy of the ground state. Besides semiconducting double 
QDs, this effect has been observed in only a few other 
systems: in metallic dots containing magnetic impurities 
and in scanning tunnelling microscopy experiments 
with two interacting magnetic atoms [12]. In addition 
to the fundamental interest in studying the Kondo effect 
on single impurities, some of those works have been moti- 
vated by the general goal of addressing and manipulating 
single spins with possible application in quantum infor- 
mation processing [l3j. Particularly promising systems 
towards this direction are magnetic molecules in which 
the exchange coupling between the spins can be manipu- 
lated by charging the molecule with an electric field [14[ . 

In this letter, the Kondo effect is used to characterize 
the exchange interaction inside a single-molecule transis- 
tor (SMT) containing a spin transition (ST) molecule. 
This type of compound supports a transition between 
a low-spin (LS) and a high-spin (HS) state, which can 
be induced by a change in temperature, pressure, mag- 
netic field or by light irradiation [15j. We have studied 
a ST molecule containing a Fe 2+ metal ion connected 
to two ligands : [Fe-(L) 2 ] 2+ [L=4'-(4" '-pyridyl)-l,2':6T'- 
bis-(pyrazolyl) pyridine] (referred to as [Fe-(L)2] 2+ ) (see 



Fig. 1(a)). This molecule shows a reversible thermally 
driven spin transition at 286 K in its bulk crystalline form 
[l6| . This type of ST molecule is particularly suited for 
studies in a SMT geometry since the charge state of the 
molecule can be modified by applying a potential on the 
gate electrode. According to recent ab-initio calculations 
|17| . this charging may drive the ST from LS to HS in 
analogy to charge transfer ST compounds [l8| . 
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FIG. 1: Schematic representation of: (a) the molecule [Fe- 
(L)2] 2+ [L=4'-(4"'-pyridyl)-l,2':6'l"-bis-(pyrazolyl)pyridine] 
and (b) the the transistor geometry. 

To connect the molecule, we fabricate nanometer sized 
gaps by a controlled electromigration process on a thin 
Au wire deposited on top of an oxidized Al gate [l9| . 
Electromigration and the subsequent self breaking of the 
wire is carried out at room temperature in a solution con- 
taining the molecule; the solvent is subsequently evapo- 
rated. If a molecule is trapped in the gap a three terminal 
device is formed in which the molecule connects via tun- 
nel barriers to source and drain electrodes and couples 
electrostatically to an external gate (see Fig. 1(b)). 

We present in Figs. 2(a) and 3(a) the differential con- 
ductance over bias and gate voltage obtained at T = 
1.5 K for two samples, A and B. The data exhibits classi- 
cal behavior of single-electron transport through a nano- 
object with a single relevant energy level. The extracted 
addition energies E ac ^ >100 meV and gate couplings 
j3 = Cg/Ctot ~ 0.02 (C g : gate-dot capacitance; Ctot- 
total capacitance of the dot) are typical of transport 
through a single molecule [19]. The tunnel couplings for 
samples A and B can be estimated from the Lorentzian 
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broadening T of the Coulomb peaks : Fa ~60 meV and 
T# ^10 meV. A ZBP, characteristic of the Kondo effect, 
is observed for both samples in a pronounced regime of 
gate voltages. For sample A, the Kondo temperature ex- 
tracted from the temperature dependence of the zero-bias 
conductance of the Kondo peak [2o| presented in Fig. 5 
is Tka ~ 70 K. For sample B, a temperature Tkb ~ 25 
K is deduced from its full width at half height. 




-4-3-2-1 1 2 3 4 




1 1.5 2 2.5 3 3.5 -10 10 

Vg(V) V(mV) 



FIG. 2: (a) Differential conductance of sample A versus 
source-drain (Vb) and gate (V^) voltages, (b) Zoom of the 
differential conductance plot, (c) Differential conductance 
traces (from top to bottom) versus Vb taken at gates voltages 
of respectively 2, 2.4, 2.9 and 3.3V. The black (red) traces 
correspond to gates voltages on the right (left) hand side of 
the diagonal line. 

In addition to this typical SMT behavior we observe 
a particular feature visible as a diagonal line (indicated 
by a red arrow) in the differential conductance plots of 
Figs. 2(a,b) and 3(a,b). In Fig. 3(a) on the right hand 
side of this line, a shift AV g of the Coulomb diamond 
edges in gate voltage is present, which finds a natural 
interpretation as a response to a modification of the elec- 
trostatic environment. This feature can be reproduced 
by a Coulomb blockade model which is based on a sce- 
nario with two parallel dots: current flows through one 
QD, the "transfer dot"' in the presence of a second dot, 
the "spectator dot" (see Fig. 4(a)). This second dot is 
electrostatically coupled to the first dot and as its tun- 
nel coupling to one of the electrodes is very weak; there 
is no net current through it. However when the gate 
voltage is increased this dot can be charged by an addi- 
tional electron, thereby, modifying the electrostatic en- 
vironment of the other dot via their mutual capacitance 
Cm- This charging induces an horizontal shift of the 
Coulomb diamond edges by AV g . Fig. 4(b) shows the re- 




FIG. 3: (a) Differential conductance of sample A versus 
source-drain (Vb) and gate (V g ) voltages, (b) Zoom of the dif- 
ferential conductance plot (c) Differential conductance traces 
(from a to e) versus Vb taken at gate voltages of respectively 
0.8, 0.9, 1.15, 1.2 and 1.25 V. The black (red) traces corre- 
sponds to gates voltages on the right (left) hand side of the 
diagonal line. 



suit of a rate equation calculation of two single-level QDs 
capacitively-coupled (see supplementary information for 
more details), which mimics the situation for sample B. 
We will discuss later the most likely physical origin of 
this double dot scenario. 

We first focus on the low-bias regime. A zoom of this 
region is shown in Fig. 2(b) and 3(b). For sample A the 
ZBP is clearly split on the right hand side of the diago- 
nal line and two peaks appear at ±eT4 = = ±2.75 
meV (see Fig. 2(c)). The temperature dependence of 
the ZBP and split ZBP are presented in Fig. 5(a,b). In 
case of the split ZBP, the zero-bias conductance G(0) 
exhibits a nonmonotonic temperature dependence (Fig. 
5(c)). For sample B the splitting is less pronounced and is 
visible as a reduction of the zero bias conductance (Fig. 
3(c)). When the gate voltage is further increased the 
splitting becomes more pronounced due to the decrease 
of the Kondo temperature as we go away from the de- 
generacy point [2l|. For the trace indicated as "e" in 
Fig. 3(c) we extract a splitting Ab ^0.5 meV for sam- 
ple B. Split ZBP as well as nonmonotonic temperature 
dependence of G(0) are reminiscent of spin-spin interac- 
tions competing with Kondo screening, as observed in 
coupled QDs systems 0, [To| and described theoretically 
in Refs. [22|, |23|. With this in mind one might expect 
that the split ZBP is induced by interactions between the 
two spins of the double dot also in our case. This expla- 
nation would be consistent with the observed transport 
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FIG. 4: (a) Schematic representation of the double dot sys- 
tem, (b) Simulation of the differential conductance ver- 
sus source-drain (V&) and gate (V^) voltages calculated for 
T = 4.2 K. (c) and (d) Schematic representation of molecule- 
electrode contacts realizing respectivly the intrinsic (Ri) and 
extrinric (R e ) double-dot scenario. 



characteristics in Fig. 2(a) and 3(a), if we assume the oc- 
cupation of the dots to be as represented schematically in 
Fig. 4(b). On the left hand side (region I), the transfer 
dot (1) has an even occupancy (as confirmed by the ab- 
sence of Kondo peak in this region in Fig. 3(a)). When 
the gate voltage is increased (region II), the transfer dot 
is charged by an additional electron, leading to an odd 
occupancy, St =1/2 and the observation of a Kondo peak 
(Fig. 2(a) and 3(a)). In region III the spectator dot is in 
turn charged leading to an odd occupancy St=S s =l/2 of 
both dots, if an even occupancy is assumed for the spec- 
tator dot on the left hand side of the diagram. A sizable 
spin interaction between the transfer and the spectator 
dot could then explain the experimentally observed peak 
splitting. However, as we will argue below the estimated 
exchange interaction between the two dots is typically 
too small in order to account for the experimentally ob- 
served value Aa,b =0.5-3 meV. Here the spin transition 
will come in. 

Two possible physical realizations of the double-dot 
scenario, and R e , are conceivable (see Fig. 4(c) and 
4(d)). In the first realization, , we make explicit 
reference to the internal molecular structure of our [Fe- 
(L)2] 2+ molecule. To support this scenario, we have car- 
ried out ab-inito density functional theory (DFT) calcu- 
lations [I?} (see also supplementary information) on the 
molecule. The two bi-pyrazolyl- pyridine ligands support 
7r-systems, which are oriented perpendicularly to each 
other by the octahedral coordination around the Fe 2+ 
metal ion. Therefore, left and right ligand states hy- 
bridize only weakly with each other. In particular, for 
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FIG. 5: (a) and (b) Temperature dependence of differen- 
tial conductance traces of sample A versus Vb taken respec- 
tively for V g =2.2V and Vg=3.2V. (c) Temperature depen- 
dence of the zero bias conductance in (a) and (b) . The 
red line represents the fit obtained from the functional form 
G{T) = G(0)[1 + (2 1/s - l)(T/T K ) 2 ] s [HHH using s = 0.22 
and Tkb=70 K. 



the situation of double reduction (i.e., when two excess 
electrons are added) our calculations confirm that each 
electron is localized on one ligand only, with very lit- 
tle excess charge penetrating the region connecting the 
two ligands near the Fe 2+ ion. In the asymmetric situ- 
ation where one ligand is connected with two electrodes 
via 7r-ligand metal-electrode half-sandwich type binding 
interactions [24J, and the other one only to a single elec- 
trode (as represented in Fig. 4(c)), the double dot system 
sketched in Fig. 4(a) is indeed realized. In the alterna- 
tive way to produce a double dot system, R e , one could 
imagine an extrinsic spectator dot, that does not form a 
chemical bond with the original molecule, e.g. a second 
molecule or a small metallic grain (see Fig. 4(d)). 

We start the discussion with the scenario R^ as 
illustrated in Fig. 4 (c). A first estimate for the 
ligand-ligand exchange coupling (with the molecule be- 
ing in the LS state) may be obtained using J^ _L ~ 
£> 2 U (£>=t/EHOMO* : relative energy splitting between 
symmetric/anti-symmetric combinations of ligand or- 
bit als; U: effective interaction energy in the near ion 
region of wavefunction overlap, roughly approximated 
by the e^-t2^-splitting: U< A oct ). We have typically 
q w0.1-l%, U^5 eV, so Ji <0.01-lmeV. This rough es- 
timate is fully supported by our ab initio calculations 
P3] (CASSCF level, see supplementary information) that 
predict for the LS state a weak ligand-ligand ferromag- 
netic exchange coupling J^ls ~0-095 meV. Essentially, 
this low value is due to the suppression of the tunneling 
matrix elements, t, between the two ligand systems. It is 
only weakly modified if a small variation of the 90° an- 
gle between the ligands is taken into account (0.15meV 
change at 20°). To further support the theoretical argu- 
ments, we can compare these values with observed ex- 
change splittings for LS transition metal complexes [HI, 
which are in the jaeV range. In conclusion, according 
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to our calculations, the scenario when assuming a LS 
state leads to an exchange coupling too small to explain 
the observed splitting Aa,b- 

Now, in the second case involving an extrinsic dot, R e , 
the direct tunneling between the spectator and the trans- 
fer dot [U [26| should be expected to be even smaller than 
in the scenario since they are not connected via a 
chemical bond. Hence, also the associated exchange cou- 
pling will be very small. Of course, in principle there is 
an additional indirect interaction due to the RKKY cou- 
pling via the leads [10|, An accurate estimation of this 
splitting is difficult to obtain as it depends strongly on the 
spatial configuration and geometry of the dots and elec- 
trodes. However, we argue that the associated exchange 
coupling should anyways be very weak. This is because 
coherent tunneling processes connecting the two dots 
have to take electrons across two barriers: from the spec- 
tator dot into the leads and then from the leads into the 
transfer dot (and the same way back to establish coher- 
ence) . The parametrical estimate (" golden rule" type) for 
the associated coupling is J^ KKY « T t T s / 'A ts / \r ts ^F) 2 
(T t ,s : electrode induced broadening of the molecular lev- 
els of the transfer and spectator dot; A st : mismatch of 
the two coupled energy levels; r ts kF' spatial distance of 
the two contact points of the two dots with the electrodes 
measured in units of the Fermi- wavelength) . Assuming 
T M wlO meV, A ts ~ 100 meV and r ts k F « 3 - 10 typ- 
ical values of J^ KKY lie between 0.01-0.1 meV. Hence 
we conclude: an explanation of the observed split ZBP, 
which would be based on a conventional spin-spin inter- 
action inside a double-dot seems hardly consistent with 
the values of the observed splittings Aa,b- 

So far, our analysis was based on the assumption that 
the molecule is in the LS state. However, a detailed the- 
oretical analysis suggests that the LS state may give way 
against a HS state when the [Fe-(L)2] 2+ molecule is dou- 
bly charged [I?}- Accordingly, the molecule could be in 
its LS state in region I and II, while it would be in the 
HS state in region III. In HS the Fe 2+ core possesses 
a spin 5 = 2 which is antiferromagnetically coupled to 
the spins on the ligands. The ground state is thus a 
spin triplet [ItJ [Hj Due to the spin-orbit interaction the 
triplet splits into a Kramer's doublet and a singlet state. 
In the optimized high-spin geometry of our molecule, we 
find that this zero field splitting, AE, is about 0.16 meV. 
This value will be larger if distortions due to incorpora- 
tion in the planar device geometry are present and our 
preliminary calculations suggests that it can reach values 
of several meV. We thus expect that among all the pos- 
sible scenarios, Rf s gives rise to the strongest splittings 
which approach the experimental values of a few meV. 

In summary we have observed a double-quantum-dot 
behavior in combination with a split Kondo peak. We 
suggest that the split ZBP indicates a charged-induced 
spin- transit ion inside the molecule. To reach a definitive 
conclusion in favor of this high-spin scenario, further ex- 



periments, which could involve, e.g., a detailed measure- 
ment of the Kondo splitting in magnetic fields, would be 
needed. The observation of such a spin-transition would 
be an important step towards the development of molec- 
ular spintronics, as it demonstrates the possibility of ad- 
dressing and controlling with an electric field individual 
spin states engineered at the molecular scale. 
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